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For the calculation of polarization transfer observables for quasielastic scat- 



tering of protons on nuclei, a formalism in the context of the Relativistic 



^% ' Plane Wave Impulse Approximation is developed, in which the interaction 

5-j ' 



matrix is expanded in terms of a complete set of 44 independent invariant 

amplitudes. A boson-exchange model is used to predict the 39 amplitudes 

which were omitted in the formerly used five-term parameterization (the SP- 

VAT form) of the nucleon-nucleon scattering matrix. Use of the complete set 

of amplitudes eliminates the arbitrariness of the five-term representation. 
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I. INTRODUCTION 

Quasielastic scattering of protons on nuclei is an attractive phenomenon for the study 
of the basic nucleon-nucleon interaction in the nuclear medium, because it exhibits the ap- 
proximate behavior of the scattering of a nucleon on only one nucleon of the target nucleus. 
Quasielastic scattering has been modeled by the Relativistic Plane Wave Impulse Approxi- 
mation [ffl which considers it as a single-step process, whereby the projectile interacts with 
only one nucleon of the target nucleus, while the rest of the nucleons remain inert. The well 
known and outstanding success of the original RPWIA was its prediction of the analyzing 
power for 40 Ca(p, p ) and 208 Pb(p, p ) at 500 MeV; a case in which all non-relativistic models 
failed [g. 

In the RPWIA approach, the description of the initial and final free particle in the 
medium is based on mean- field theory, as described by Serot and Walecka in Ref. ||. In 
the RPWIA model the associated Dirac plane waves have their free nucleon mass decreased 
by the real part of the average nuclear scalar field to yield an effective nucleon mass. The 
values of the effective masses serve as an indicator of the nuclear medium effects on the NN 
interaction. 

In former theoretical studies of scattering [[|||[1| the nucleon-nucleon scattering matrix 
(F) was parameterized in terms of the five Fermi covariants which is commonly referred to as 
the SPVAT form of F or the IAl model. It should be stressed, however, that eventhough the 
SPVAT form gave reasonable results for elastic and quasileastic scattering observables, it is 
in principle not correct, since as was first pointed out in Ref. ||, a five-term representation of 
the relativistic NN scattering matrix is necessarily ambiguous. In addition, Tjon and Wallace 
|7j have shown that a general Lorentz invariant representation of F (referred to as the IA2 
model) contains additional terms that cannot be neglected. The IA2 representation of F 
contains, in fact, 44 independent invariant amplitudes, instead of the previously used five, 
which are consistent with parity and time-reversal invariance as well as charge symmetry 
together with the on-mass-shell condition for the external nucleons. Comparison to the 



limited data available with subsequent and more refined calculations P,|8|-[T0| have also 
revealed that quasielastic (p,p ) and (p, ft) prefer different five-term representations of F: 
The (p, ft) data favor a pseudovector irNN coupling, whereas the (p, p ) data are consistent 
with a pseudoscalar term for the nNN vertex. Therefore, the most basic question which 
has to be addressed is the representation of the NN scattering matrix. 

In the current application of the RPWIA to quasielastic scattering, the following com- 
ponents play a key role: 

1. The amplitudes in the basic two-nucleon interaction, which are partly determined from 
free NN scattering data and partly from a solution of the Bethe-Salpeter equation 
employing a meson-exchange model for the NN force. 

2. The Lorentz covariant set constructed from the Dirac matrices which serves as a rep- 
resentation for F. 

3. The effective nucleon mass for both projectile and target nucleons interacting in the 
nuclear medium. 

In this paper a theoretical formalism is presented for the calculation of calculate polariza- 
tion transfer observables for quasielastic proton-nucleus scattering using a general Lorentz 
invariant representation of the NN scattering matrix [] IJ; a systematic survey of the predic- 
tive power of the model compared to data will be presented in a future paper. By adhering to 
the symplifying features of the RPWIA, one can focus on the basic NN interaction without 
introducing additional complications. A complete and unambiguous expansion of F allows 
for a correct incorporation of medium effects (within the RPWIA framework) and there- 
fore one can distinguish between experimental results which are genuine medium effects and 
those which arise due to other effects not taken into account by the RPWIA. In Section |Tl 
we briefly review the RPWIA and also discuss the ambiguities of the SPVAT form of F. 
In Section [111] the general Lorentz invariant representation of F is discussed. Section |IV] 
presents the transformation from invariant amplitudes to effective amplitudes, while in Sec- 



tion [V] expressions for the spin observables are derived in terms of the effective amplitudes. 
A calculation of complete sets of spin observables, based on the IA2 model, for quasielastic 
40 Ca(p, p ) scattering at 500 MeV is presented in Section [VT]. Section |V11| summarizes the 
main aspects of this paper. 

II. RELATIVISTIC PLANE WAVE IMPULSE APPROXIMATION 

Complete sets of spin observables [] (P, A y , Dm, D s / S , D s n, D[/ s , D nn ) for quasielastic 
(p, p ) and (p, ft) scattering are calculated within a relativistic framework using the Rela- 
tivistic Plane Wave Impulse Approximation (RPWIA) |l[. The RPWIA models quasielastic 
scattering as a single-step process whereby the projectile knocks out a single bound nucleon 
from the nucleus. The rest of the nucleons are assumed to remain inert, but their effect is 
taken into account in that the free mass of the projectile and target nucleons are shifted 
to effective masses, Mi and M 2 respectively. In the context of the Walecka model || the 
effective masses can both be calculated microscopically as follows: For the projectile, 

M 1 = M+< S(f ) > 

where M is the free nucleon mass, < S(f ) > is the average value over the whole nucleus 
of the real part of the scalar potential found by weighting it with T(b), the transmission 
probability through the nucleus at an impact parameter b, and with the nuclear density p(r) 
for the specific nucleus ||. The effective mass of the target nucleon is determined from 

M 2 = M - g s < 4>(f ) > 

where g s is the scalar meson coupling constant and 4>(r ) is the scalar field for the specific 
nucleus with the averaging done as described above. Values of Mi and M 2 for specific nuclei 
and incident laboratory energies can be found in Table II of Ref. [H]. Experimental data 



1 The spin observables are defined in Section [V] 



seem to suggest that the spin observables are target independent |T2|,|T3j, and therefore we 
assume, as a first step, a Fermi-gas approximation for the target nucleus. The RPWIA 
therefore reduces quasielastic proton-nucleus scattering to a two-body scattering process 
with Dirac spinors (containing effective nucleon masses) describing the external nucleons. A 
graphical representation of the scattering process is depicted in Fig. [TJ. 
Refering to Fig. [I], the projectile Dirac spinor is given by: 



U(p 1 ,M 1 ,s l 



E{ + M 1 



2Mi 



[*i 



o-'Pl 



(2.1) 



where E{ = pi + Mf and the spinor is normalized to U(pi, M 1 , s) U(pi, M ly s') = 5 ss i. 
Similar expressions exist for the other three spinors labelled by p 2 , k\ and k 2 . The following 
four-momenta are also defined: 



where pf = Mf and kf 
spin projection operator, 



p\ = (El ft), p* 2 = (E* 2 ,p 2 ), 

k\ = [E*;,k x ), k* 2 = (E*',k 2 ) 

- Mf (i = 1, 2). For handling the polarization, one requires the 



P(n) = -(I 2 + a-n) 



(2.2) 



for the direction n where I 2 is the 2x2 unit matrix. In the basis of Pauli spinors, <fi(n) for 
spin direction n, we have 



p(fi) = cj)(n)(j)\n). 



(2.3) 



Defining 



U{p 1 ,M 1 ,s i ) = U j >{p 1 ,M 1 ,s i ) 1 D 



where the convention of Ref. |14| is used for the gamma-matrices, the Lorentz invariant 
matrix element for the scattering process depicted in Fig. p] is given by: 



M = [U(h, M u s f ) ® U(k 2 , M 2 , 4)] F p/(pi, M l5 s< ) ® [/(p 2 , M 2 , s 2 )] (2.4) 

where F is the 16 x 16 nucleon-nucleon scattering matrix. The question arises as to what 
form of F is to be used in Eq. Q2.4J) assuming parity and time-reversal invariance as well 
as charge symmetry. Once a choice of F has been made, analytical expressions for all 
spin observables, namely the unpolarized double differential cross section, the analyzing 



power and the polarization transfer observables can obtained from Eq. ( |2.4j ). All previous 
calculations |H,|15|,|8|-|i0f of spin observables for quasielastic proton-nucleus scattering have 
been parameterized F in terms of only the five Fermi covariants: 

F = F s {h <8> h) + F P (j 5 <g> 7 5 ) + F v {^ <g> 7m ) + 

F A ( 7 V ® 7 5 7m) + M^ V ® ^) (2-5) 

where the latter is commonly called the SPVAT form or IAl representation of F. The 



amplitudes, F L (L = S, P, V, A, T) are obtained by fitting to free NN scattering data [|15 
This procedure, however, does not uniquely fix the form of the matrix F . To see this we 
note that the pseudoscalar covariant, PS = 7 5 <S> 7 5 , has exactly the same matrix elements 
between positive energy free mass Dirac spinors (M x = M 2 = M) as the pseudovector 
covariant, PV = ^- <g) ^-, i.e., 

'U X {M) ® U 2 (M)] [PV - PS] [U ± (M) ® Z7 2 (M)] = 0. 

This is called the equivalence theorem [0. We can therefore replace PS with PV in Eq. 



( p.5| ) without altering the amplitudes, .F&. Eventhough these two representations are equiv- 
alent on-shell (p 2 = M 2 ) they will give different results when sandwiched between positive 
energy Dirac spinors containing an effective nucleon mass, since then matrix elements be- 
tween negative energy states now also enter. This is because the effective mass spinor can 
always be expanded in a free mass basis: 

U(pi, Mi, Si) = auU(pi, M, Si ) + ayV(pi, M, s t ) 



where V is the negative energy Dirac spinor [14|. There also exists the relation 



MiM 2 , , , N 

Mpv = ~J^T M PS (2-6) 

where M.ps and M.py are the contribution of the pseudoscalar covariant and pseudovector 
covariant respectively to the invariant matrix element given by Eq. ( p.4| ). Note that in 



Eq. Q2.6|), the pseudovector covariant is PV = %rh~ ® qm~i but where q = p\ — k{ 



2M ^ 1M 

k 2 — p 2 , i.e., the momenta are on-mass-shell with respect to the effective masses, Mi and 
M 2 . In the equivalence theorem, the momenta must be on- mass-shell with respect to the free 
mass. The above equality has been used in Refs. J8"HT0f to investigate the sensitivity of the 
spin observables to the difference between using a pseudoscalar covariant or a pseudovector 
covariant. The ambiguity which is inherent in any five-term or incomplete representation of 
F (such as the IAl representation) was first pointed out in Ref. |J. Tjon and Wallace have 
developed a general Lorentz invariant representation of F. The formalism can be found in 
Refs. |7|,|TTf and is applied to elastic proton- nucleus scattering in Refs. p"?l-|T9|j. We will refer 



to this as the IA2 representation of F and discuss, in the next section, its application to 
quasielastic proton-nucleus scattering. 

III. IA2 REPRESENTATION OF F APPLIED TO QUASIELASTIC 
PROTON-NUCLEUS SCATTERING 



From Eqs (3.1) and (3.18) in Ref. |19| the IA2 representation of F is given by: 

13 

P = E E F Z lP ' lP2P ' 2 t A Pi (kil M) ® A p , (k 2 ; M)} K n 
p\p' x \pip' 2 n=1 
[A pi (p 1 ;M)®A P2 (p 2 ;M)} (3.1) 

where M refers to the free nucleon mass. Henceforth the notation 

{p} = PlP'l,P2P' 2 

will be used. In Eq. ( |3.1| ), F^ (n — 1 — 13) are the invariant amplitudes for each rho-spin 
sector (which is defined by the rho-spin labels, pip[; p 2 p 2 , where p = ± ), A p (p,M) are 
covariant projection operators given by: 



AMM) = elJ^L 

p(E 7 ° -p-7) + M 

2M l ' ' 

where E 2 = p 2 + M 2 , and K n (n = 1 — 13) are kinematic covariants constructed from the 

Dirac matrices: 

K\ = S = I 4 <g> I 4 

K 2 = P = 7 5 ® 7 5 
^3 = V = 7 M ® 7/, 
iT 4 = A = 7 V ® 7 5 7m 
iT 5 = T = ^ (8) <v 
K 6 = QnAh®in 
K 7 = Q22A1 11 ® h) 
^8 = Qu lA ,(7 6 ®7V) 

^9 = Q22 lAt (7 5 ®7V) 

^io = Qi2, M (/ 4 ®7 M )£ 
^n=Q 2 i, M (7 M ®/4)5 

^i2 = Qi2, At (7 5 ®7V)5 



K i3 = Q 2 iAr'f®r)s 



where 



Qv> = - M 3 M with j^ = fe x and pf, = h- 

With each combination of rho-spin labels { p\p'\Pip'<i } is associated a pair (ij) to index a 
specific rho-spin sector (or subclass): See Table I of ]nj. For example {+ + ++} = (11) and 

{H 1 — } = (22). Parity and time-reversal invariance, together with charge symmetry and 

the on-mass-shell condition for external nucleons, lead to F being completely specified by 
44 independent invariant amplitudes fill . Five amplitudes in subclass F 11 are completely 

8 



specified by fitting to physical free NN scattering data and are therefore identical to the 
SPVAT amplitudes in the IAl representation of F. The remaining 39 off-shell amplitudes 
(contained in subclasses F 12 to F 44 ) are obtained by solving the Bethe-Salpeter equation in 
a three-dimensional quasi-potential reduction |20| , |2~l"|j , with pure pseudovector pion-nucleon 
coupling, to determine a complete set of helicity amplitudes. The invariant amplitudes are 
related via matrix equations to the helicity amplitudes fTTJ. The IA2 representation is a 



complete and unambiguous expansion of F since covariants cannot be added or changed 
arbitrarily without violating the above-mentioned symmetries. Amplitudes which are solely 
determined by physical scattering data are isolated in subclass F n while the remaining 
amplitudes are determined by solving a dynamical equation, the Bethe-Salpeter equation 
using a meson-exchange model for the NN force. 

From Eq. ( p.l|) four cases concerning the combination of projectile and target nucleon 



masses can be distinguished: 

1. No medium effect (Mi = M 2 = M): In this case only subclass F 11 will contribute to 
the invariant scattering amplitude. 

2. Projectile Relativity (Mi 7^ M; M2 = M): Contributions to the invariant scattering 
amplitude arise from F ll ,F 21 ,F 31 and F 41 where the latter three subclasses require 
at least projectile relativity for a contribution. 

3. Target Relativity (Mi = M; M 2 7^ M): Contributions to the invariant scattering 
amplitude arise from F ll ,F l2 ,F 13 and F 14 where the latter three subclasses require 
at least target relativity for a contribution. 

4. Target and projectile relativity (Mi 7^ M; M 2 7^ M): Now all subclasses will contribute 
to the invariant scattering amplitude but F 22 , F 23 , F 24 , F 32 , F 33 , F 34 , F 42 , F 43 and 
F 44 require at least projectile and target relativity for a contribution. 



From Eq. (|3.1|) we see that medium effects can never occur in subclass F 11 due to the ac- 



companying positive energy projection operators. Medium effects in the IA2 representation 

9 



of F arise only due to off-shell amplitudes (which are contained in the subclasses F 12 to 
F 44 ). This is in contrast to IAl where medium effects are included only in subclass F . 
One can now substitute Eq. (|3.1|) into Eq. Q2.4J) and proceed from there to calculate the 
spin observables in terms of | M. \ 2 which is directly related to the invariant amplitudes F^ . 
We will, however, not follow this direct approach due to the following reasons: 



1. Following the standard procedure (see Ref. |14[ for example) one finds that \M\ con 



tains traces over at least eight gamma matrices. The number of gamma matrices 
increase as the covariants become more complicated. Since the number of terms gen- 

erated by such a trace is given by — ^ (where N refers to the number of gamma 

matrices), and since there is a double sum over the rho-spin sectors, a very large 
number of terms will occur. 

2. Since we are applying a relativistic formalism to a Nuclear Physics problem, it might be 
more instructive to rewrite the NN scattering matrix in a form which is more familiar 
to traditional Nuclear Physics. 



We will therefore follow a similar approach as in Ref. |22] where an effective t-matrix is 
derived which is a 4 x 4 matrix, but which still contains all the information coming from the 
relativistic analysis. From Eq. ( |2.1| ) we can write: 



I7(pi, Mi, Si ) = (j±) 2 « + (pi, Mi)0(*) 

where, as a 4 x 2 matrix: 



u+(pi,M 1 
Similarly 






2E{ 



'2 



(3.3) 



where, as a 2 x 4 matrix: 
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u + (pi,Mi)=« +t (pi,Mi)7°. 
u p (p) (where p — ±) contains no reference to the spin and is normalized to 

u p '\p'p)u p (pp) = 5 plp . 
In terms of u + the invariant matrix element, (Eq. |2.4| ) is given by 

M = [ ^ip J ^(sf)u + (h, M x ) ® 0t( 8 ;) fl +(fe, M 2 )] F 
[0(sO« + (pi, MO <g> 0(s 2 )« + (p 2 , M 2 )]. 

Use of the identity AC ® BD = (A <g> B)(C ® D), where <S> refers to the usual Kronecker 
product, leads to the expression 

M = ( EI m?mP ) ^ s /) ® <^)] [^ + (^i> MO ® «+(4 M a )]F 

[u + (pi,M0 ®w + (p 2 ,M 2 )][0(sO ® 0(s 2 )]- (3.4) 

Defining the effective t-matrix as: 

t = [u + (k h MO ® w + (fe, M 2 )]F[w+(pi, MO <8> w + G» 2 , M 2 )] (3.5) 

and ^ = [(F*F 2 *F/F 2 *V^iA#, Eq. O becomes: 

M = 01 [0 f O/) ® <f>\s' 2 )]t[<f>(si) ® 0(«2)]- (3.6) 

Since £ is a 4 x 4 matrix it can be expanded in terms of a basis constructed from the Pauli 
matrices and the momenta of the scattering process. Define the three-momentum transfer 
q = pi — k\ = k 2 — f>2, the average momentum p a = \{pi + k\) and a vector orthogonal 
to both q and p a , N = q x p a = p 1 x k\. Note that 

q-p a = ~{Pi - h 2 ). 

For quasielastic scattering \p\\ ^ \k\\ and therefore q and p a are not orthogonal, however, 
N ■ q = N ■ p a =0. Assuming only parity invariance, t can be written in terms of a set of 
eight linearly independent matrices in the spin of the two interacting nucleons: 

11 



v (1) - ■ T 

Xl — 1 2, 


(2) 

Xi 


X? = N-a, 


(2) 

x\ -- 


X ( 3 ] = N-a, 


(2) 
X3 


v (1) -- r 

Xi — 1 2, 


(2) 

xl -- 


(1) -, -, 

xl = q ■ <*, 


(2) 

xh 


(i) _ -> 

Xe = Pa' V, 


(2) 


(1) - -, 
Xi = q ■ a, 


(2) 

Xi ~- 



i=T.b n {x^®X^) (3.7) 

71=1 

where 

"" = h- 

N-a. 

= h- 

N-a. 

q-a. 

Pa -o- 

Pa ■&■ 

(1) -. -, (2) -> -. 

X& = Pa-a, xl = q ■ o. 
In the next section the invariant amplitudes, Ffi (n = 1 — 13) are transformed to a set of 
eight effective amplitudes b n (n — 1 — 8), and expressions for the spin observables are derived 
in terms of the effective amplitudes. 

IV. TRANSFORMATION FROM THE INVARIANT AMPLITUDES TO THE 

EFFECTIVE AMPLITUDES. 

Expressions for the effective amplitudes of t are now derived: Taking the trace of Eq. 
(0) yields 

6 1 = i'&[t]. (4.1) 

Multiply Eq. fl3.7| ) with (N ■ a <g) N ■ a) and take the trace of the resulting equation. Since 
N ■ q = N ■ p" a = there will be no contribution from the last four terms of Eq. (|3.7|) and 
therefore 

1 



4N 
where 



MN,N) (4.2) 



Y 2 (a, b) = Tr[(o • a ® b ■ a)i 
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Similar arguments lead to 

1 



AN 
and 

1 



-Tr[(N ■ a ® I 2 )i] (4.3) 



■Tr[(I 2 ®N-o)i]. (4.4) 



AN 2 

Following the same reasoning as above, one can also derive a set of four coupled equations 
relating the amplitudes 65, be, b-j and 6s- A set of coupled equations arise since the vectors 
q and p a are not orthogonal for quasielastic scattering, (i.e. q-p a 7^ ). The solutions are: 



r 2 Y 2 (p a ,p a ) - r 2 r 3 Y 2 (p a ,q) - r 2 r 3 Y 2 (q,p a ) + r 2 Y 2 (q,q 

M 2 (rir 3 — r 2 y 
riY 2 (p a ,p a ) - nr 2 Y 2 (p a ,q) - rir 2 Y 2 (q,p a ) + r 2 2 Y 2 (q,q 



, '2 1 2\Pa,Pa) - i2'3*2KFa, V ) ~ > 2> ?^ 2\H, Pa) T '3 12\H,H ) ,. -.-. 

°s tt^ttz zi^ ^- b > 



6 ~ M 2 ( ri r 3 - r 2 ) 2 [ ' 



-nr 2 Y 2 (p a ,p a ) + r 2 2 Y 2 (p a ,q) + rir 3 Y 2 (q,p a ) - r 2 r 3 Y 2 (q,q) 

M 2 (r 1 r 3 - r 2 y 
-rir 2 Y 2 (p a ,p a ) + rir 3 Y 2 (p a ,q) + r 2 2 Y 2 (q,p a ) - r 2 r 3 Y 2 (q,q) 

M 2 (r 1 r 3 - r 2 y 



— 1 \> 2 L 2 \y a ,ya) T 1 2 i 2 yp a ,q J T il'3*2\<d,Pa) '2'3*2\q,q ) , . „, 

b 7 = " I757Z _ Zl^ l 4 -0 

2 ) 

, 1 .1'::!' 2WV- I', 1 I - ' U :<,'- 2\.l><i-<l 1 - r 2 Y 2 yq,p a , 1 - 2 i y, 1 ;M>/.'/ ; r 

"8 = ^ToTZ Z Z~2Y2 ^-°> 



where 



2q 



n = 


M ' 






2p a - 


q 


•2 


M 




r 3 = 


IP' 





and 



M 

The next step is to derive an expression for the t-matrix which is convenient for use in the 
calculation of the traces which determine the effective amplitudes. Substitution of Eq. ( p.l|) 
into Eq. (|3.5|) leads to 

i = Y,lt F n } [ f P ; (h, M, MO ® f p , 2 (k 2 , M, M 2 )] K n 

[T P1 (pi, M, MO ® T P2 (p 2 , M, M 2 )] (4.9) 

where we have introduced the 4x2 T-matrices defined as: 
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T p {p, M, M*) = A p (p, M) u + (p, M*). 



(4.10) 



In Eq. ( [4.1 0| ) M* denotes an effective mass and Eq. ( |3.3| ) has been generalized to: 

( h<P(pi,M*) ) 



u + (p,M*) 
Eq. ( |OTD reduces to Eq. (^3J) if we set 



a 



PX(P,M*) 



(4.11) 



' E*(p x ) + Mi 
, 2E*(pi) , 



where E*(p 1 ) = Jpf + Mf and 



x(pi,Mi) = 0(pi,Mi) (£*(pi) + Mi))" 1 . 



We can obtain an explicit expression for T p as follows: From Eq. (|4.11|) we can write 

\ 



u + (p,M* 



(l, 
\ 



2 <P(p } M 




+ 



/ 







a-px(p,M*) 



V 



J 



and therefore 



u + (p, M*) = (/>(p, M*) (ei ® I 2 ) + X (P, M*) (e 2 ® o* • p ) 



(4.12) 



where 



ei 



M 



vV 



and e 2 



M 



V 1 / 



with ejej 



that 



5y. To write the p-spin projection operator in (2 x 2) <g> (2 x 2) form we recall 



'/, 0^ 



7 



P-l 



-J 2 



cr 3 ® J 2 and 



p • a 
-p ■ a 



\ 



2<r 2 <E> p ■ a. 



(4.13) 



(4.14) 
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Substitution of Eqs. (gl5|) and ( gig ) into Eq. Q) leads to 



^-^(0-3(81/2) —(o 2 ®p-o) + -I 

2M V ; 2M V y 2 



A p (p,M) = ^(a 3 ®/ 2 ) - -il^®^) + -{h®h). (4.15) 



Substitution of Eqs. Q4.12J ) and ( |4. 15| ) into Eq. ( |4.10| ), and using the properties of the Pauli 

matrices, allows one to write: 



r p (p,M,M*) = J2hf(p,M,M*) [e,(g)A(p)] (4.16) 



i=i 



where 



MP 




with 



hPtf, M, M*) = £§0(p, M*) - ^- X (P, M*) + h(p, M*) and 



2ikT v "' ' 2M /vyri ' 2 



h p 2 \p,M,M*) = j^mM*) - §fX(p,M*) + lx(p,M*). 
Similar steps lead to 

T p (p, M, M*) = J2 3?{P, M, AT) \e\ <g> A(p)} (4.17) 



8=1 



where 

jjP (p, M,M*) = h p 1] {p, M, M* ) and 

jf (p,M,M*) = -hf\p,M,M*). 

To calculate the effective amplitudes, the contribution of each covariant to the trace relations 
must be determined, t is calculated from Eq. ( [4.9| ) using the explicit forms of Y p . (p, M, M*) 
and r pi (p,M,M*) in Eqs. ( gl§ and (pTQ . Use is then made of Eqs. flO]) to Q to 
determine the effective amplitudes, b\ to b$- This procedure requires the calculation of traces 
of a set of matrices £,, where i = 1 — 46. With each covariant is associated a set of t- matrices 
of which the trace must be taken. We wrote a program in the MATHEMATICA computer 

15 



language to do the required trace algebra. The effective amplitudes are linear functions of 



i#», i.e. 



bi = bi({F^ } }) where i = 1-8. (4.18) 

The isospin zero (isospin one) effective amplitudes are obtained by substituting the isospin 
zero (isospin one) invariant amplitudes into Eq. ( 14.181) . 



V. EXPRESSIONS FOR SPIN OBSERVABLES IN TERMS OF THE EFFECTIVE 

AMPLITUDES 

In this section expressions are derived for the unpolarized cross section, the analyzing 
power and the polarization transfer observables in terms of the effective amplitudes b n for 
both (p, p ) and (p, ft) scattering. 

Working in the nucleon-nucleon laboratory frame, the spin in the incident beam direction 
is described in terms of three orthogonal unit vectors (l,s,h), where I is along the beam 
direction, s lies perpendicular and to the side of I in the scattering plane, and the normal 
unit vector is ft = I X s. Similarly the spin of the final beam is described in terms of 
(l',s\h). 

1. Unpolarized double differential cross section 

For the scattering process in Fig. [I] one can write down the expression for the differ- 



ential cross section, da [14 



\v x -v 2 \\ E{ El ) \E\ E* 2 J ^ X 2 ' 2 

d 3 ki d 3 k 2 , , j. 2 



(2vr) 3 (2tt) 



-\M\< 



As we consider the quasielastic scattering at energies much higher than the interaction 
energies amongst the target nucleons, we assume the latter to be practically non- 
interacting. Therefore, the momentum distribution of these nucleons can be obtained 
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in a Fermi-gas model. Following the same arguments as in Ref. [|]] allows one to write 
down the following expression for the double differential cross section: 

' d\p 2 \ d<f> \p2\f(pl,P2,h f )\M\ 2 



where 



and 



dQ' 1 dE[ \q\E{ h 



f(Pl,P*2, k f) 



2/1 #2 



M{M£ 



167r 3 £;3[(^. p *_ M 2 M 2)2]i 



Pr, 



q lu* 
2 ~T 



4M| 



In Eq. ( |5.2|) , q^ is the four-momentum transfer g M = (a;*, q ) where 



(5.1) 



(5.2) 



uj* = E{ — E{ = E 2 — Ei and q = Pi — k\ = k 2 — f>2- 

To obtain the Fermi momentum kf, the required effective density is calculated in an 
eikonal approximation as shown in Ref. |IJ. More refined values of kf for specific target 
nuclei can be found in Table II of Ref. ||. Define the function 

r"(p 1 ,p 2 ,k h k 2 ) = "£ £ \m\ 2 . 



*2- 



Substitution of Eq. ( |3.7D into Eq. ( |3.6| ) leads to 



M = gi'£b n tfisAxZWs 



•Hs'2)xf^S 2 ) 



n=X 



and therefore one can write 



T"(p 1 ,p 2 ,k u k 2 )=g 2 1 V b* m b n TT[ X ^X^}Tr[x^X^} 



m,n=l 



An explicit expression for T (pi,p 2 ,ki,k 2 ) is given in the appendix. To obtain the 
unpolarized cross section, one sums over the initial spin and average over the final spin 
which leads to: 
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d\p 2 \ d(j) \p2\f[Pi,P 2 ,kf) 



dQ[dE[ J unpol A\q\E{ Jpmin 

T"(pi,---k 2 ). (5.3) 



For (p,p ) scattering 



da \ (p,p ' ] IhlE' r k f 

1 d|p 2 | c?0 |P2|/(p*,P2,fc/) 



d^dE'J unpol 4|g|^'^ 

(z B// r"(pi, • ■ .4 {&,(/ = 1)}) + iv e// r"(p 1; • • • k 2 , {br}) 



where 



i(&,(J = 0)+ 6,(1 = 1)). 



For the charge-exchange reaction (p, n) 



da \ M IhlE', fh 

d\p 2 \ d(f> \P2\f(p 1 ,P 2 ,kf) 



dQ[dE[ J unpol A\q\E\ Jp mi 

eff 



N eff r"(p 1 ,---k 2 ,{b?- ex }) 



where the charge-exchange amplitudes are defined as: 

6f— = ^[6,(7 = 1) -6,(1 = 0)]. 

The quantities Z e ff and N e ff are defined in Ref. [l|] and values for specific targets are 
given in Table II of Ref. g. 

2. Analyzing power 

The definition of the analyzing power is given in terms of polarized double differential 
cross sections as: 

where, for example, 

da I fci| Ei /•% » ^fV-M 

V ,rn/ ( g /) = |^| / / d\p 2 \ d(j) |p 2 |/(Pl,P 25 fc /)o r ( S /) 



dQ 1 dE l \q\E* Jp min " 2 



is averaged over incident spin directions s$, and the target particles' initial and final 
spin as contained in the factor: 

f'(p 1 ,...,4%)= E l-^l 2 

= nx^P{s f )x^Ax^X^] (5.5) 



where use was made of Eqs. (|2.2| ) and Q2.3|) . A calculation of the traces in Eq. ( |5.5|) 
shows that Y (pi, ...,s/) has the following structure: 

f'(pi,...,s/) = /i(pi, ...,&) + f2(Pi,-,k 2 )N- s f + 

f 3 (p h ...,k 2 )Pa-(qxs f ). (5.6) 

Defining the combination function 

r'(pi,...,fc 2 ,s/) = f'(pi,...,fc 2 ,s/) - f'(pi,„.,£ 2 , -S/). (5.7) 

and using Eq. ( |5.6|) in Eq. (|5.7|) yields 

r '(s/) = 2(/ 2 (pi, ..., k 2 )N ■ s f + f 3 (p u ...,k 2 )p a - (qx s f ). 

The explicit forms of the functions / 2 and f 3 can be inferred from Eq. (|A2|) in the 
appendix. If Sf = n then 

N ■ h = pxki sin^ 

and 

Pa- (q* h) = -pikismd L . 

The analyzing power for the (p, p ) reaction is given by 

A r vx = /pj; < „d|^l#|P2l/(pi,..,P^fc/)(Z e// r / (n,{b i (J=l)}) + jv e// r'(n,{6ri)) 
P ' P " J^L dlftl #1 ft|/(*/) (^e//r"({6,(/ = 1)}) + N eff T"({br})) 

and the analyzing power for the (p, n) reaction is given by 

= pjjg|jf|, (5 . 8) 

/ P l^|p2|#|p 2 |/(A ;/ )r({6f- ea '}) 
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Since a (p,n) reaction implies that the incident proton could only have scattered off 
a neutron, we set Z e ff = and therefore N e ff appears as a common factor in the 
numerator and denominator and cancels out, which means that N e ff does not appear 
in Eq. 



3. Polarization transfer observables 

The polarization transfer observables are defined in terms of linear combinations of 
polarized double differential cross sections as follows: 



,/ 7T -,/ 



A ' j ;sfe(*> h) + ikw(~^ h) + Jfefr, -St) + ?ksrX-ii> -*/) ' (5 ' 9) 



In Eq. ( |5.9| ) a typical polarized differential cross section is: 

jo' j f' (*"*/) = 7^7777 / °^ 2 1^ ^/(Pi^&fkf^s/) ( 5 - 10 ) 

a\l 1 ahj 1 \q\E* J p mi n * 

where 

f(p 1 ,...,k 2 ,s i ,s f ) = J2 \ M \ 2 

S2,S 2 

= 9l E ^fcn[Tr(P(^)xWP(l/)x«)][Tr(^) X ^)] 

m,n=l 

= /i(pi, ..., fc 2 ) + -4*1 • s< + A 2 ■ s f + 

(k ■ A 3 )(s f ■ A A ) + (Si ■ s f )(A 6 ■ At) + h ■ (s f x A 5 ) 

with Ai functions of only the three-momenta, p\ to hi of which the explicit form can 



be inferred from Eq. (AT). Define again, now dictated by the form of Eq. ( |5.9|) , a 
function: 

r(pi, ..., k 2 , Si, sf) = r(si, §f) - f (-Bi, sf) - f(si, -§f) + 
r(-§i,-s f ). 

= 4[(s< ■ A 3 )(s f ■ A 4 ) + Si ■ (sf x A 5 ) + 

(si-s f )(A 6 -A 7 )}. (5.11) 
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The explicit expression for T contains various kinematical parameters which are pre- 
sented in the first column of Table |. The other columns contain the values of these 
quantities in the laboratory frame, for each polarization transfer observable: 9 refers 
to the laboratory scattering angle, pi = \p\\ and k\ = \ki\. Use of Eqs. ( |5.10| ), ( |5.11| ) 
and (]5.3|) lead to 



I p k i m d\p 2 \ dcj> \p 2 \f(pi, P * 2 ,k f ) r (a, ...,&) 



1 J rkf 11 -» I J / I -* I /•/ * * J \ tV/ / -» I \ 



The polarization transfer observables for the (p, p ) reaction are given by 

D ^ = Ipj m d\P2\ d4 \p2\f(kf)^Z eff T(s u s f ,{b t (I= 1)}) + 4^^,^, {6?™})) 

*'* " SiL dm d<t> lP2\f(kf) {ZeffT"({k(I = 1)}) + AT e// r"({&r e })) 

and the corrsponding observables for the (p, n) reaction are given by 



Once again, as in Eq. ( |5.8|) , the effective number of neutrons does not appear in Eq. 

(EH- 



Although the primary aim of this paper is to present the theoretical formalism for calcu- 
lating quasielastic proton-nucleus polarization transfer observables, in the next section we 
give a brief glimpse of the predictive power of the formalism by applying it to quasielastic 
40 Ca(p, p ) scattering at 500 MeV A systematic study of the predictive power of the model, 
as well as a comparison to IAl-based predictions, will be presented in a future paper. 

VI. RESULTS 

Before presenting the results we mention the numerical checks that were performed to 
verify that the transformation from invariant amplitudes F^, to effective amplitudes b n 
was carried out correctly and that the expressions for the spin observables in terms of the 
effective amplitudes are indeed correct. For Mi = M 2 = M only subclass F 11 contributes 
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to the invariant matrix element and the IA2 representation is therefore equivalent to the 
SPVAT form of F, We therefore verified that our expressions for the spin observables in 
terms of the effective amplitudes give exactly the same numerical result as the corresponding 
expressions in Ref. [jXJ] which contain only the five SPVAT amplitudes. This confirms that 
the transformation to effective amplitudes has been carried out correctly for only the SPVAT 
covariants. To verify the transformation for covariants K§ to K13, we derived expressions for 
the spin observables directly for each individual covariant K§ to K\^. This involves traces 
over Dirac matrices (as opposed to the trace algebra involving Pauli matrices presented in 
this paper) and provides a non-trivial check for the transformation involving covariants K§ 
to Ki%. The fact that two independent ways give numerically the same result for all spin 
observables confirms the correctness of the transformation to effective amplitudes and the 
expressions for the spin observables derived in this paper. 

The formalism in the previous sections is now applied to quasielastic 40 Ca(jo, p ) scatter- 
ing at an incident laboratory kinetic energy of 500 MeV and a laboratory scattering angle of 
19°. In the original calculation of Horowitz and Murdock in Ref. |TJ, it was found that the 
use of an effective mass for both the projectile and target nucleons moved the theoretical 



calculation closer to the data [23| and below the free mass calculation for A y . This was 



referred to as the quenching effect in the analyzing power and claimed to be a "relativistic 
signature" . In Ref. [l[] the SPVAT parameterization of F was used. Fig Q shows the results 
employing the IA2 representation of F. The solid line represents the calculation using an 
effective mass for the projectile and target nucleons with ^- = 0.892 and ^ = 0.817 
taken from Table II in Ref. M for 40 Ca at T/ a & = 500 MeV. The dashed line is the free 



mass calculation. The data are from Ref. [23|. We notice that the quenching effect in A y is 



very small compared to Fig. 6 of Ref. [I] over the entire energy range. The result is that 
the IA2 calculation does not describe A y as well as the IAl calculation of Ref. [JJ. For the 
other observables, the effective mass and the free mass calculation do equally well. This is 
in contrast to the result in Ref. |T| where the Zfyj's only preferred a free mass calculation. 
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VII. SUMMARY 

We have presented a theoretical formalism to calculate polarization transfer observables 
for quasielastic proton-nucleus scattering using a general Lorentz invariant representation 
of the nucleon-nucleon scattering matrix. In this way we avoid the ambiguities which are 
inherent in the previously-used five-term representation (the SPVAT form) of F . In the 
process we have derived an effective t-matrix, which is a 4 x 4 matrix and therefore more 
familiar to Nuclear Physics, but which still contains all the information coming from the 
relativistic analysis. This necessitates the transformation from the 44 invariant amplitudes 
to a set of eight effective amplitudes as well as the derivation of new expressions for the 
spin observables in terms of the effective amplitudes. By staying within the framework of 
the Relativistic Plane Wave Impulse Approximation (with its many simplifying features) | 
and using a general Lorentz invariant representation of F, allows us to do an unambiguous 
investigation of medium effects via quasielastic proton-nucleus scattering. The first applica- 
tion of the formalism to the reaction 40 Ca(p, p ) at I} & = 500 MeV and 9i a b = 19° shows 
that the IA2 representation of F does not lead to such strong medium effects in any of the 
spin observables, in contrast to the results in Ref. [1[] where the medium effect was most 
noticable in A y . There it was also found that the use of an effective mass for the projectile 
and target nucleons led to the theoretical calculation being closer to the data than the free 
mass calculation. The IA2 representation is consistent with data, however, in that it predicts 
little medium effect in any of the spin observables, eventhough the prediction of A y is now 
a little poorer than before. In a subsequent paper, a systematic study of spin observables, 
using the IA2 representation of F, will be presented for both quasielastic (p, p ) and (p, n) 
data. 



2 all of which are motivated by experimental data on the spin observables 
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APPENDIX: EXPLICIT EXPRESSIONS FOR SPIN OBSERVABLES IN TERMS 
OF EFFECTIVE AMPLITUDES a, L (i = 1-8). 

In this appendix we present explicit expressions for the quantities T , T and T in terms 
of the effective amplitudes a; L which are related as follows to the effective amplitudes bo 

61 = di 

, a 2 

bi = — -aj for i = 3,4 and 

m 2 



1 
bi = —di for 2 = 5,6, 7, 



m? 



where m denotes the free nucleon mass. 



— r (pi,p2,ki,k 2 ) 
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'4Im(a 2 ) 2 4Rc(a 2 ) 2 



4Im(ai) z + 4Rc(ai) z + (N ■ N)' 

\ III" III" I 

A Al ( Um(a 3 ) 2 4Rc( fl3 ) 2 4Im(a 4 ) 2 4Rc( fl4 ) 2 ^ , ( AJm(a 6 ) 2 4Rc(a 6 ) 2 ^ 2 

N ■ N [ + + + + + ( Pa ■ Pa ) 

8Im(a 6 )Im(a7) 8Rc(ae) Rc(a7) 8Im(o6) Imfag) 8Rc(a 6 ) Re(a 8 )\ „ _ 

: 1 : : 1- -\ : I Pa ' Pa Pa ' Q + 



mr rrr m^ m 



1 






1111(05) Imfffle) 8Rc(a5) Rc(ag) 8Im(o7)Im(o8) 8Re(a7) Rc(a 8 )\ , _ _,2 

; 1 : ■ 1- : ■ h : I (j> a ■ Q ) 



to" 1 nv* nv* m"* ) 

8Im(a 5 )Im(a7) 8Rc(a 5 ) Rc(a 7 ) 8Im(a 5 ) Im(ag) 8Re(a 5 ) Rc(a 8 )\ 



4 , Pa-qq-q + 



' 4Im( ffl5 ) 2 4Rc(a 5 ) 2 ^ ^ 2 ( 4lm(a 7 ) 2 4Rc(a 7 ) 2 , 4Im(a 8 ) 2 4Rc(a 8 ) 2 ^ _ .. „ _ 
3 + 3 \ KQ-9) + 3 + 3 + 3 + 3 Pa-p a q-q (Al) 
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-^r (pi,...,fc 2 ,s/) 



m 



-4Re(a 2 )Im(a 4 )JV-iV 4Im(a 2 ) Re(a 4 ) TV • TV 4Re(ai)Im(o 3 ) 4Im(oi) Re(a 3 )\ 



N-S f + 



4Re(a 6 )Im(a 7 )p a ■ p a 4Im(a 6 ) Re(a 7 )p a ■ p a 4Re(a 5 ) Im(a 6 )p a • g 4Im(a 5 ) Re(a e )p a ■ q 



Pa- (<fx Sf) 

m m m m 

4Rc(a 7 ) Im(a 8 )p a • Q 4Im(a7) Re(a 8 )p a • q 4Rc(as) Im(a 8 ) q • q 4Im(a5) Re(as) q ■ q 



(A2) 



-^T(p 1 ,...,k2,s i ,Sf) = 

f 2lm(a 2 ) 2 N ■ N 2Rc(a 2 ) 2 N ■ N 2Im(a 3 ) 2 2Rc(a 3 ) 2 \ - - „ f 2lm(a 2 )Im(a 4 ) N ■ N 

s 1 5 1 Z 1 Z N ■ Sf N ■ Si+ \ - h 

\ in" in" m 4 m 4 J \ rn° 

2Re(a 2 )Rc(a 4 )/V-/V 2Im(ai) Im(a 3 ) 2Re(oi) Re(a 3 )\ -t A . _ , /2Im(a 6 ) 2 p a • p a 

5 5 2 N -{ Sl X Sf)+p a ■ S fPa ■ S, 

m° to z m z I \ m 4 

2Rc(a 6 ) 2 p a -pa 4Im(a 6 )Im(a 8 )p a ■ g 4Rc(a 6 ) Rc(a 8 )p a ■ q 2Im(a 8 ) 2 q- q 2Rc(a 8 ) 2 q- q \ 
to 4 to 4 to 4 to 4 ' to 4 J 

_ „ /2Im(a 6 )Im(a 7 )p a ■ p a 2Re(a 6 ) Re(a 7 )p a • p a 2Im(a 5 ) Im(a 6 )p a • g 2Re(a 5 ) Re(a 6 )p a ■ q 

Pa 'Si + j + j + j + 

\ m m 4 m 4 to 4 

2Im(a 7 ) Im(a 8 )p a • g 2Re(a7) Rc(a 8 )p a • q 2Im(a5) Im(a 8 ) q- q 2Rc(as) Re(a 8 ) q ■ q\ _ „ 

4 1 1 1 1 1 4 1 ■ s f + 

to 4 m 4 m 4 to 4 J 

_ „ f2lm(a 6 )lm(a 7 )p a -p a 2Re(a 6 ) Re(a 7 )p a • p a 2Im(a 5 ) lm(a e )p a ■ q 2Re(a 5 )Re(a e )p a ■ q 

Pa-Sf 1 + 4 + J + 1 + 

\ to 4 m m m 4 

2Im(a7) Im(a 8 )p a • q 2Rc(a 7 ) Re(a 8 )p a • q 2Im(a5) Im(a 8 ) q ■ q 2Re(a5) Re(as) q- q 



q- Si + 
in- in- in- j 

2Im(a 7 ) 2 p a ■ p a 2Rc(a 7 ) 2 p a ■ p a 4Im(a 5 )Im(a 7 )p a -g 4Rc(a 5 ) Re(a 7 )p a • q 2Im(a 5 ) 2 g-g 
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4 4 4 

mr rrr mr 
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FIGURES 
FIG. 1. Two-body scattering process with momentum, mass and spin labels for the external 

nucleons. F is the 16 x 16 nucleon-nucleon scattering matrix in the two-nucleon spin space, pi and 

ki (i = 1 2) represent the four-momenta of the particles respectively. M\ and Mi denote the 

effective masses of the projectile and target nucleons respectively. Si, S2, s/ and s 2 are the spin 

four-vectors for each particle. 

FIG. 2. Spin observables for a range of transferred energy uo over the quasielastic peak for 
inclusive proton scattering from 40 Ca at 500 MeV and 9^ = 19°. The centroid of the quasielastic 
peak is at u; ~ 63 MeV. Data are from Ref. [23]. The solid represents the IA2 calculation and the 
dashed line represents the free mass calculation. 
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TABLES 
TABLE I. Expressions for kinematical quantities containing Sj and/or §f for each non-zero 
polarization transfer observable. 



kinematical 
quantity 


Dm 


D s i s 


J-^nn 


D sn 


D Vs 


q- Si 


P\ — k\ cos 9 


—k\ sin 9 





P\ — k\ cos 9 


—k\ sin 9 


q-s f 


P\ cos 9 — k\ 


—pi sin 9 





—p\ sin 9 


p\ cos 9 — k\ 


Pa 'Si 


^ {pi + k\ cos 9) 


2&i sin 9 





\{pi + ki cos 9) 


■yk\ sin 6* 


Pa'Sf 


\{p\ cos 9 + fci) 


— \p\ sin9 





— \p\ sin9 


\{j>\ cos 9 + /ci) 


N-Si 








P\k\ sin 9 








N-s f 








P\k\ sin 9 








Si-Sf 


cos 9 


cos# 


1 


— sin# 


sin# 


N ■ (sj x sf) 


P\k\ sin 9 


p\k\s\n 9 





P\k\ cos 9 sin 9 


—p\k\ cos 9 sin 
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